Spin flip probability of electron in a uniform magnetic field 
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The probability of a spin flip of an electron is calculated. It is assumed that the electron resides 
in a uniform magnetic field and interacts with an incoming electromagnetic pulse. The scattering 
matrix is constructed and the time needed to flip the spin is calculated. 

PACS numbers: 42.50.Ct 
Keywords: spin flip 



INTRODUCTION 

Controlling the spin state of an electron is becoming 
an increasingly important issue. Spintronics continues 
to advance and managing the spin of the electrons is a 
vital concern. Quantum information seeks long lived cu- 
bit states that allow reversible transfer and control. [l[ , 
Moreover, in generalized theories of gravity the spin gives 
rise to a new field, called torsion, which may someday be 
measured from spin flipping techniques. [3], 4] 

A review of experiments involving spin flips in solids 
has been given, [5j and more recently experimental in- 
vestigations of nonlinear optical effects for spin polarized 
electrons in semiconductors has been undertaken. Q Spin 
flip by a dc voltage in a quantum dot was also reported. 

For these reasons it is important to know the effect 
of an electromagnetic wave on the spin state. A classi- 
cal investigation was undertaken by Gover,Q who stud- 
ied electron beams interacting with an electromagnetic 
pulse. While the main subject was spin resonant radia- 
tive emission, it was noted the spin flip is enhanced when 
coherent light is used. A related problem dates back to 
Brown, [9] who studied the geonium "atom," which is a 
single particle placed in a Penning trap. In this case the 
magnetic field is not uniform. 

In this Letter we compute the quantum mechanical 
probability of an electron spin flip induced by an elec- 
tromagnetic pulse. In particular we consider an electron 
trapped in a uniform magnetic field, but otherwise free. 
Now we assume an electromagnetic pulse is incident on 
the electron and ask, what is the spin flip probability. 

We start with the transition amplitude, defined by 



Sfi = -ie I d^xip p° Aaijji 



(1) 



where / and i denote final and initial states, A a is the 
electromagnetic potential, and for now h = 1 = c. 

Normally the final and initial states are free particle 
wave functions, but that is not appropriate here. For 



one thing, from a philosophical view, it is not sensible to 
describe about the spin of an electron without a fiducial 
field by which one can measure the spin. In order to 
quantify the spin, there must a field, usually a magnetic 
field. Thus we consider the scattering of a bound state to 
another bound state, so the asymptotic wave functions 
should be those of an electron in a magnetic field, which 
is assumed to be constant. 

The wave function for an electron trapped in a uniform 
magnetic field in the z direction was found long ago and 
is given byfioj 

^„ = C f e- iiEnt - p >- p ' z) e-^/ 2 u n (2) 
where the momentum terms are eigenvalues, 
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where the energy eigenvalues are 
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where h n = N n H n , H n are the Hermite polynomials, 
N n = 1/ ■\/2"n!y / 7r, and by definition the Hermite poly- 
nomial with a negative subscript is zero. The u n are 
given below. 

One precautionary note is given before the calculation 
is presented. The cross section for spin flip is inherently 
small because the electromagnetic interaction mixes the 
"large" part of the Dirac spinor with the "small part" 
(whereas, for example, a chiral term mixes large with 
large) . Thus other terms not involving spin flip are much 
larger, in general. These terms include momentum trans- 
fer effects, where the electron acquires a velocity (as in 
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Compton scattering in lowest order). Thus, an electron 
at rest trapped in a uniform magnetic field can, besides 
spinflipping, acquire a net velocity from the interaction. 
However, if the electron is bound to a solid, then the elec- 
tron will be held in place. However, in this case the the 
interaction Hamiltonian will differ than that used here, 
in that there is an additional potential to be included. 

In order to get an exact spin flip probability one must 
use the exact bound state wave functions and potential 
of a bound electron in an external field. However, such 
conditions are not known exactly, in general, and so the 
calculation used here may be taken as an estimate, pro- 
vided the internal fields are small compared to the exter- 
nal field. For example, in a typical paramagnetic salt the 
local magnetic field of the order of 50 mT, which is small 
compared to the 4 T field used here. 

Now we will calculate the probability of a spin flip 
starting with the scattering matrix 



Sfi = lim t 



< ipf\4>i > 



(6) 



As mentioned above, we are considering the scattering 
between two (quasi) bound states, but the analysis fol- 
lows along the usual lines. We start with the Dirac equa- 
tion for a particle bound in a uniform magnetic field, with 
potential A a , that experiences an incoming electromag- 
netic wave with potential A a , 



(ij' 7 d (T -e^(A a + A cr ))^ = mil). (7) 
The propagator is defined according to 

(z7 ff 9 a - ej* A a - m)A = 5(x - y) (8) 
so that we have the general solution 

4>(x) = <f>(x) - i [ d 4 yAV(y)^(y) (9) 



where V = ej a A a and cf> is the solution to the homo- 
geneous equation. A full treatment of the propagator is 
not needed here since we are only interested in the first 
order effects, so that 
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The essential difference between this and ([T]) is that in 
(fTUf the wave functions correspond to a free particle but 
in the above the functions are those of an electron in a 
magnetic field. 

Now we consider the spin flip such that the final state 
is spin up and the initial state is spin down. In this case 
we have, 



'n-1 





p (h n ^ /(E f +m) 
\ -y/2neBh n /(E f + to) J 



(11) 



and 
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In the above all of the h n are h n (£), where £ is the 
dimensionless argument given above (£ = U eB j hey — 
—p x c/ ' \eBhc in cgs). 

It is assumed that we have box normalization, but only 
in a two dimensional box. This is so because ip is a func- 
tion of y. Thus we may impose periodic boundary con- 
ditions in x and z, but not in y. In essence, in order to 
obtain the Hermite polynomials, it is assumed the solu- 
tion must vanish at infinity. Imposing different boundary 
conditions would yield solutions that differ than those 
given above. 

The sides are given by L x and L z . It may be noted that 
the gauge freedom in the choice of the electromagnetic 
potential translates to a freedom in the choice of y or x, 
or a combination of the two. 

To begin the calculation, let us write (fTUl) as 



Sfi 
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M fi e-^ Eit - p >- p >^ (13) 
where Lq = u)/ eE and the matrix element is 



M fi = uJ 7 VA, 



(14) 



which becomes 
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Our main interest lies in spin nipping and not on trans- 
lation. Therefore we focus here on the terms that are 
independent of the initial and final momentum (other- 
wise, if these eigenvalues are zero, we see that the matrix 
element vanishes). For this reason we take the electro- 
magnetic potential to be of the form A a = {0, 0, 0, ^.3}. 

We see from our choice of the potential, the matrix 
elements occur with Hermite polynomials with different 
indices. Due to the orthogonality we see that this term 
would vanish unless the potential is a function y. For 
these reasons, assuming the electromagnetic field is a 
Gaussian pulse, we take 



E = Eq cos(ky — wt)e 



(16) 



where in the following we drop the subscript on the am- 
plitude and d is a dimensionless constant that essentially 
determines the number of wavelengths in the pulse. We 
assume that d » 1 so that the potential may be written 
as 



A 3 « -l sin ( k y - ^)e-(^) 2 EE ^f, (17) 
CO LU 

which gives the electric field in the x direction, the mag- 
netic field in the — z direction, and the wave propagating 
in the y direction. 

Using this in (|13|) we can perform the x and z integrals 
yielding delta functions, 



Sfi 



Lo 



(18) 



where X = 2nS(pl — p x ) and Y = 2-k8{p{ — p\) and 



a = J dtdye-i 2 fe^f-^Uf^Asu, 



(19) 



Now, we define a± — a + — a_ so that, after integrating 
over time, 



a± 



2iuj 



K 



dye 



\Fle~B 



n i — n f 



E; + m 



Ef + m 



(20) 



where K = -(dA+/2w) 2 , A ± = E f - E t =F w, and b = 

ikL^l + A+fu). 

For below, we define L x = 1/y/eB and 4> — Px/VeB- 
To prepare for the final spatial integration we assume a 
photon is absorbed, so we arc interested in a + which is 
simply called a from here on, and write this as, 



a = N 2 L l J d£e- e+btL Mf 



(21) 



where 



> — dc 
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(22) 



We make the substitution /j, = £ — 6/2, complete the 
square, and use the identity, 



H n (x + y) = p(j^H k (x)(2 y y 



— k 



(23) 



where 



are the binomial coefficients, so that 



a = N 3 ^EE^ 

n k f )(lty h k(0K(0b n f- k b n *- k -R 



(24) 



where N3 = \/2eBhcN2Lie b / 4 and R is defined as the 
entire term preceding it with rif and n, replaced by 
and nj_i. 

The integrals may be performed without trouble since 
the h n are orthonormal with the weight function e - ^ . 
The result is 



L (E f 
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where 
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Now we define P as 



P 



\Sfi\ 2 dp 



(27) 



where p is density of states. Earlier we adopted two di- 
mensional box normalization in the x and z directions, 
so the phase space of this volume is 



p = dp x dp z 



{2-Kf 



(28) 



Performing the phase space integral we finally have, 



where Q = b 2 /A - d 2 (A+) 2 /4o 
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If P is less than one then we may interpret this as the 
probability that the electron suffers a spin flip. However, 
as defined P may be greater than one, which physically 
corresponds to the fact that the electron will eventually 
flip back, and then again, and so on. In this case P 
represents the number of times the particle flips. In many 
cases we are interested in investigating the parameters for 
which P = 1. For example, let us calculate how long it 
takes for the particle to flip. 

Suppose we consider the transition rii — 0, n,f — 1. Let 
us further assume that the electromagnetic field is tuned 
to the transition energy, and take p x = 0. With this the 
result reduces to P = ird 2 E 2 /B 2 , which may written in 
terms of the intensity of the electromagnetic wave, / as 
(in cgs) 



8^ 
cB 2 



(30) 



Consider the question, how long must the pulse last 
in order to flip the spin. Earlier we defined d and the 
number of wavelengths in our pulse, so that dX = ct. 
Now we solve for the time, t = T, such that P = 1. The 
result is 



T = 



One should note this result holds at resonance only, off 
resonance the scattering matrix decreases exponentially 
with time, and the corresponding spin flip time increases 
quickly. 

This result seems to be compatible with a result of Hu 
et. al., who determine the spin flip time of an electron 
to be greater than 0.1 ^s. To emulate those results we 
consider a magnetic field of 4.0 T in a field of 200 mW. 
To estimate the intensity we assume a distance of 30 cm 
and a standard 16 dB gain antenna, which gives T = 
27 fis. However, in that experiment, they considered a 
low density of electrons CdTe in quantum wells, so the 
comparison cannot be given too much weight. 



Alternatively one may ask what intensity is required 
to flip an electron in a given time interval. Using the 
parameters given above, for example, it would take about 
5 W cm -2 to flip the electron in one microsecond. 

In the above it was assumed that d » 1, which means 
there are many wavelengths in the pulse. It is natural to 
ask what the result for small d is. To investigate the 
effects of a short pulse let us take the potential to be 
the Gaussian form of width d. This corresponds to a 
field that looks like a single wavelength of the order of 
d in length. The calculation may be carried through as 
above, and the result is that the time needed to flip a 
pulse is reduced by a factor of 2d 4 . 

In summary, the probability of the spin flip of an elec- 
tron in a uniform magnetic field was calculated in the 
case that an electromagnetic pulse was incident on the 
electron. At resonance the time it takes to flip the spin 
was estimated. 
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